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ABSTRACT 


Imagine a signal acquisition system composed of a 
number of receivers or sensors concurrently scanning the 
same domain for signals. It is reasonable to expect that 
different signals will each be detected by a different 
subset of receivers over the scanning period. Using the 
data collected from the receivers, it is possible to es- 
timate the total signal population size including those 
Signals not detected by any receiver. Additionally, it 
is possible to estimate the probability each individual 
receiver detects signals. Several estimators are devel- 
oped for these quantities in the context of a model 
designed to represent the signal detection process. This 
model forms the basis for a simulation conducted to analyze 


the behavior of the estimators over a variety of conditions. 
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ir LOWE TION 


Imagine transmission devices emitting signals in the 
electromagnetic spectrum at unknown times and of unknown 
quantity. Further, imagine there are receivers concur- 
rently scanning the spectrum in an effort to detect the 
transmissions. Suppose there are k receivers in the sys- 
tem and a record is kept of each signal detected by each 
receiver. It is reasonable to expect that some signals 
will be detected by only one receiver, some by more than 
one receiver, and some signals will not be detected at 
all. Assume each signal is identifiable by its charac- 
teristics making it possible to determine whether detec- 
tions by several receivers were actually detections of a 
common signal. The purpose of this thesis is to demon- 
Strate how data collected by the acquisition system can 
be used to estimate the total number of signals exposed 
to the system including the number of signals not de- 
tected by any receiver. It would, additionally, be 
interesting to know the probability each receiver in the 
system detects signals, and this problem will also be 
discussed. 

It would stand to reason that the detection process 
is some function of the receiver's capability to detect a 


Signal as well as the signal's ability to be detected. 





One such functional relationship (see Barr [Ref. 1] ) will 
bemproposed in the context of a model developed to repre- 
sent eee ional detection process. Various estimation 
principles such as maximum likelihood, method of moments, 
and least squares will be used to develop estimators for 
the total signal population size and the receiver intercept 
probabilities. 

A development of several estimators will be shown in 
detail and their applicability to this problem discussed. 
These estimation procedures will then be incorporated into 
a general model constructed to represent the detection 
process. This model supports the conduct of a computer 
Simulation which also includes the generation of signals 
with varying signal strength and simulates detections of 
the signals by the receiver system. The resulting esti- 
mates are compared to the actual number of signals gener- 
ated and the expected receiver intercept probabilities to 
provide some insight on the quality of the estimators. 

The estimators are exposed to various conditions of the 


Signal process and their behavior analyzed. 





II. GENERAL MODEL 


In a statistical approach to modelling the physical 
phenomenon of signal detection, it is desirable to develop 
a model of sufficient generality such that the factors 
associated with radio and radiowave theory can be broadly 
parameterized. This generality has a distinct advantage in 
that it is not necessary to precisely define what type of 
sensor is being used nor what necessarily constitutes a 
Signal. A detection is presumed to be some interception 
of a signal by a sensor which enables the signal to be 
identified. It may be that a signal, once detected, is 
lost and then redetected. This model similarly presumes 
that a reoccurence of a signal may also be identified. It 
is also desirable to prescribe a uniform method for counting 
Or indexing detections in a manner which is conveniently 


applied to mathematical formulae and computer code. 


ee SLGNAL DETECTION 

It is possible to postulate functional values for the 
probabilities of signal detection See be defined to 
be the probability that receiver i detects signal j, inde- 
pendent of other signals. Each receiver or sensor, regard- 
less of similarities to one another, will be peculiar in 


its relative ability to detect a given signal. Each signal 





Was) CelLtain characteristics such as intensity or fre- 
quency which affect its ability to be detected. (Depending 
on the pre eesncas . a Signal may have a variety of 
identifiable characteristics. For the purpose of this re- 
port, these characteristics are considered to be aggregated 
and called the "strength" of the signal.) Let the random 
variable V be defined as some measure of these character- 
istics where 0<V<1. Let T. (Vv) represent a specific 
function that receiver 1is "able" to detect a random sig- 
nal. Thus, a general expression, ce = eG relates 

the ability of receiver 1 to detect a signal, to the ability 
of signal j to be detected. For any future random signal 
ot, = T, (Vv). 


The form of the function considered in this model is 


I. = ii, (Yi) ania 


i ; °* V, where 1. denotes the proportion of 


Signals of strength V that receiver i would detect. It is 
easily seen that if all signals were certainly detectable, 
V would take on the value one and T, (v) = 1.. Sane) aes 
1f all receivers were of equal specification and certainly 
able to detect a given signal for any v, 1, = 1 and 

T. (v) =Secoueaiion. “Various conditions on 1. and distrib- 


mMemons Lor V will be considered. 


me METHOD OF INDEXING 
Consider that over some specified time, s (unknown) 


Signals will be scanned by k receivers in a signal intercept 
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system. For any given signal, define the random variable 2Z 
to be a k-component vector whose components GOnsitst sol 
= O if receiver i did not detect the signal and Zz. = 1 
1f receiver 1 did detect the signal. Thus, Zs is a Bernoulli 
random variable. Let S denote the sample space of Z, so § 
contains a k-dimensional vectors of zeros and ones. This 
naturally includes the zero vector which contains the un- 
observable outcome that no receiver detected the signal. 
Over some specified amount of time, for each zeS, let 
n, denote the number of times that outcome Z = z occurs and 
let n be a 2*_aimensional vector whose components are the 
n,'s. For example, inak = 2 receiver system, S = (0,0), 
(1,0), (0,1), (1,1)} . The first element, the zero vector 
is an unobserved signal. The second element represents a 
Signal detected by receiver ji only, the third by receiver 2 
only, and the last represents a signal detected by both 
receivers. Any given detected signal will appear as an 


outcome on Zé€S. Over many signals, is the number 


No=(1,0) 


of times that receiver 1 only detected the signals. The 
sum over all n, will be the total number of signals detected 


by the acquisition system. Because n is the number 


P= O00) 
of signals not detected, ae n,=s, where s was the total 
unknown number of signals present. Let N be the random 
vector with outcome n. It follows that N is distributed 


according to the multinomial law, NrwM ,(SrP), where p 
2 


eZ 





is the 2*_aimensional vector denoting the probabilities 


that a type z outcome occurs. 


Ms 


3 


Dei | ool TON PROCEDURES 


Several procedures are considered for estimating the 
total number of signals which pass through an acquisition 
system. Estimators for receiver intercept probabilities 
are developed for each approach. The maximum likelihood 
approach uses the properties of the multinomial distribu- 
tion in conjunction with numerical methods for maximizing 
the likelihood function. A method of moments procedure is 
developed which finds a solution vector to a set of simul- 
taneous equations based on the binomial characteristics of 
each receiver. Finally, a least squares approach, again 
with numerical methods, minimizes the square error between 
actual and expected detection outcomes for find estimates 
of signal population size and receiver intercept 


probabilities. 


A. MAXIMUM LIKELIHOOD/CONDITIONAL MAXIMUM LIKELIHOOD 

The concept behind this procedure is one of estimating 
the size of a multinomial population with incomplete obser- 
vations. This technique follows an approach discussed by 
Sematnanan (Ref. 2]. It is convenient, in the statistical 
sense, to imagine a set of cells representing the aes 


elements in the sample space S described in section II. A 


detection by some subset of sensors in the k-receiver system 
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Beas a count Of one to the appropriate cell. Recall that 
after many signals, the total counts in each cell will com- 
prise the 2*_dimensional vector N. N is distributed 


marecinomial, 


pa Es 


Suppose the ordering within the vector N is such that the 
number of signals which missed detection is placed cell ne 


and has probability Po: Let r denote the number of remain- 


k re 


maeecells, r = 2-1. Define t = meen. , ens me = Ss—t. 
i 2 O 
The observation of (N4,.--,0,) yields the likelihood func- 
ion L, 
n n 
1 E Sac 
oe ! 
L(s,p) s! (p,) -+e (DL) (Do) ; 


a An acai 
N,+.--np,:(s-t): 
L(s,p) may also be written as the product 


L(s,p) = L, (s,p)L, (p) 


where 
ic Ss-t 
Late, >) & = (1=p_) wp.) 
{ Sst. - a 
ae ne 
() L. (p) = t: Fy) -++ (qd) 
Nyeeeen: 
with 


a Saws)? slat are ee cee 


L, is the likelihood based on the probability of t and hence 
L. is the likelihood based on the conditional probability 

eit (Ny,eee ey) given t. The following lemma is known (see 
eee chapman (1951)). 


is, 





LEMMA: For any given p, $ = CAGeDe)) (greatest integer 
(t/1-B,)) maximizes L,(s,py), where L,(s,D)) is 
defined above. If ace = t/S' for some integer 
S', then 8 and 8-1 both maximize L,(s,p,). Other- 


wise, § is the unique maximum. 


. r 
Since aa oe 


uniquely specify Po: Hence, to find an estimate of s, it 


= 1, the sum over Das eee reels Wa 


is only necessary to maximize the conditional likelihood 
mmc t1On L.(p). More conveniently, the maximum of the 
log-likelihood function is found. Equation (1) may be 
transformed as follows: 

(2) log (L, (p)) eho gk n,log(q,) +...t nlog(q_), 
mmere K = te/(n4e...n!). 


Equation (2) can be rewritten as 


log(L,(p)) = log(k) + L& 


>_k 


i=1 


- m,log( 7; ) 7 


(t-m;) log (1- 71; ) - 


(cee (ee THD), 


1=1 

where ms denotes the sum of cell counts in S for which 
receiver i had detections (1.e. receiver i only detected, 
receiver i and receiver j both detected, etc.) and t-m, 
represents the number of non-detects by receiver i. Thus, 
it is possible to solve directly for estimates of the 
receiver intercept probabilities, 1, without solving for 
the ae components of p. 

Since K does not depend on p, it is not necessary to 


include log(K) when maximizing (2). Given that each value 
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of ns i40 is observable, it 1S possible using a numerical 
method such as the generalized reduced gradient to maximize 
the function for optimal values of TT, - The value for Po 


is clearly By = TT -f and from the lemma, $ = Cee) 


B. METHOD OF MOMENTS 
As previously defined, m denotes the total number of 
detections made by receiver i. Since each detection can be 


considered a Bernoulli trial, the random variable M; is 


distributed Binomial (s,7,). For a k-receiver system 
ees 
n, =f, 
eg 
mM. = nao 
m = $f, 


The parameters s and 1, are all unknown, hence k + 1 un- 
knowns exist in the above k equations. To obtain ak + ee 
equation, consider all the signals which were detected bv 
two or more receivers in the system. Define 
* pe ‘ 
m = . .(number of signals detected 
ies 
Simultaneously by re- 


ceivers i and j), 


If independence between receivers is assumed 


ky) _ 
E(m*) = 5s a 


1 
=- - . Mm 


me 
S 1<j 


ee 
The estimators for s and ™, thus become See mem) 7m 
a * 
and fq. mm / 2-5 25 Mm, M. - 
Wa 





ess ASlt SQUARES 
Define A to be the event that one or more receivers 

detect a given Signal. Hence, P(A) is the probability that 
the acquisition system will successfully intercept any given 
Signal. Now, m is defined to be the number of events for 
which receiver i had detections, and t, the total number of 
detections by the system. Following a procedure developed 
by Knorr [Ref. 3] , it is possible to postulate an estimator 
mor mM. as follows: 


()) ae = (m, /t) ¢ P(A) 


In this expression, m,/t denotes the proportion of signals 
detected which were observed by receiver 1. 

feeener2z CO use tne mumber of recelvers which detect a 
Signal, 


= k = = 
Za es er O02 sk 


A value for z, can be assigned for each signal intercepted. 
Further, let 13 represent the number of occasions where 
Z))= j, and let BP. Be Whe Mass “Euction for 2. It fol- 
lows that P. (j) will be the probability that exactly j re- 
ceivers will detect a given signal. A possible estimator 
for | is 


(4) p 


Nece that for j = 0, 


“~ 


(19 /¢) ¢ P(A) 


~~, 

\_j. 

ee” 
{ 


ae a 


WD 
—- 

= 

it 


cs 


Both estimators for 1; and P., (j) rely on estimating 
the value of P(A), where 0S P(A)sS 1. By setting a value to 
P(A), say P(A), and given values for ms, TJ and t, one can 
find p, (j) as a function of Pieveeiicing (4). Uciuneg the 
same mee for P(A), one can similarly compute , using (3). 
MeniS possible to find a similar value for Ee (3), denoted 
as Cee SveMaki ng use=Oof 1ts.functional ane ienentis EO 


fe - Knorr suggests the following algorithm to compute 
Pp 


a* a> iss ouster ore 


k 
~ ae k 
3. Calculate Oe (3) = ae Eon, ay 


Now form the square error of P(A), where 


2 is | aay Sa en 2 
(5) E°(P(A)) = 2F_y(P, (3) - BL G))~. 


By varying P(A) over its range, it iS possible to find a 
“least squares" estimate, O=P(A)=<1, which minimizes the 
above expression. Using DLA), Tt. can be found using (3). 
The total signal population size passing through the 


on 


acquisition system is estimated by $ = t/B(A). 


ie 





ivy, Cee eStMUEALTION 


Three methods have been discussed which provide esti- 
mators for determining receiver intercept probabilities. 
Similarly, estimators are found for the total number of 
Signals passing through an acquisition system given only a 
portion of the signals were detected. In this section, the 
techniques used to test each method will be discussed. 

The three methods of analysis were used as the basis for 
a model developed to support a computer simulation. The 
model was designed to perform two major functions necessary 
for the simulation. The first function is to generate sig- 
nals simulating the transmission of signals from some source. 
The second function of the model is to statistically relate 
each signal to each receiver to determine actual detections. 
During the simulation, each receiver must either detect or 
not detect a given signal. Thus, an accounting procedure 
was developed relating each receiver to each signal. Fi- 
nally, the model applies each method of analysis to the data 
generated and estimated values of intercept probabilities, 
ee and total signals, S, are comrared to the E (47), and 
the number of signals specified to the signal generator. 

The simulation is performed in two parts. The first 
consists of one replication of each method in each of 72 


different situations. A situation or environment is 
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comprised by varying 
fon iemrOral number of signals generated. 
2. The number of receivers in the acquisition system. 


3. The parameters specifying the relative "ability" 
of a receiver to detect a signal. 


4. A parameter specifying the relative "abilitv" of a 
Signal to be detected. 


The second part of the simulation replicates twelve of the 
above environments 100 times each, thus enabling an analysis 


of mean absolute error, bias, and variance of the estimators. 


Pee LNPUT: “MAIN ROUTINE” 

The main routine reads all required data and subsequently 
calls each subroutine in sequence. Input requirements for 
the simulation include the random number generator seed, the 
number of signals to be generated, the number of receivers 
in the system, each receiver's ability to detect a signal 
of strength V, the parameter specifying the "ability" of a 
Signal to be detected, and the method of analysis to be 
used. 

Various loops in the main routine can be established to 
test one or more methods on many data sets (environments) 


or to replicate over one data set. 


B. SIGNAL GENERATION: "SUBROUTINE GNRATE" 
Recall that T, (V) represents the specific function that 
receiver iis able to detect a random signal and that the 


Memeo: the function used (see II. A.) in this model is 
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(6) eee. = V 
In the simulation, 1. 1s a user input denoting elements of 
the probability vector specifying the ability of receiver i 
to detect any given signal. Now, suppose V is a Beta(«,1) 
random variate denoting the ability of a signal to be de- 
tected. The Beta distribution was selected due to its ease 
and simplicity of generation and its desirable behavior as 
alpha is allowed to vary. All outcomes on V are on the 
range (0,1) and can be envisioned, under certain conditions, 
to be the "probability" that a random signal v is "detect- 
able". Furthermore, it is possible to imagine a case where 
V is degenerate at 1 denoting that signal v is certainly 
detectable. Clearly, when « = 1, v becomes a uniform(0,1) 
outcome. This implies that the "ability" of signal j to 
be detected is random with an equal likelihood of being 
undetectable as detectable. Thus, by varying the « param- 
eter of a Beta(x,1) random variable, it is possible to 
specify to what degree each signal is expected to be "de- 
tectable". 

For integer «=1 and @= 1, the Beta distribution has 
the following cdf: 

Pisano, 1) = v™, Oa 7 = 1 

Simee all values of v~ lie on the interval (0,1), it is a 
Simple matter to generate random Beta («,1) variates using 


ene a 
the probability integral transform of V ~U. Let U be a 


Mecorm(0,1) random variable and va an outcome on U. 


Ao 





Then, an 


co — ol 


and ' V yl 


Because the lim ee EoneaesUertcientiy large , the 
-~ 


Beta (x,1) distribution can, in general, be used to repre- 
sent the case where TT. = 1... An « value of 100 was 
selected to simulate a signal which is virtually certain 
to be detectable. An intermediate case of x = 10 is also 
tested. 

Each signal Wa is therefore defined to be an outcome on 
er In the physical sense, « = 1 represents a system where 
all receivers, in some way, act alike or show dependency. 
For example, if an outcome on V~U(0,1) is a low value, 
nearly all the receivers in the system, regardless of their 
ability to detect the signal, will be restricted from 
doing so simply because the signal is not readily detect- 
able. Any outcome on V has this affect on the entire sys- 
tem, thus causing the receivers to show interdependence. 

On the other hand, if V is degenerate at 1 (simulated by a 
large alpha value), the detection of the signal is based 
solely on the characteristics of the receiver. In this 
manner, the receivers tend to act independently of one 
another. Hence, changes in alpha, in a restricted sense, 
give some indication of correlation through receivers. The 
case where « = 10 demonstrates an intermediate degree of 


dependence between receivers. 


Ze 





In this model, each signal v is scanned by each receiver. 
As previously mentioned, a function of the model iS 
statistically relate each receiver to each signal to deter- 
mine actual detections. T; (v) 1s computed for each receiver 
using (6) and represents the probability that receiver 1 
detects signal j}. To determine if the detection occurred, 
T, (v) is compared against a uniform(0,1) random number, 
where if 

Uieor |) T, = detection occurs 
UikOwot,) T, > no detection occurs. 

Recall that Z is a k-component vector of zeros and ones 
denoting which receivers did or did not detect a given 
Signal. Hence, Z = z2 is an outcome on the k comparisons 
above, and becomes the element in the sample space S which 
specifies which combination of receivers detected the 


sugnal. 


Sec OUNTING: ‘SUBROUTINE COUNT” 

The estimators developed in section III rely on certain 
"counts" of events related in various fashions to the re- 
ceiver system. Subroutine COUNT takes the data created by 
the signal generator and performs the necessary accounting. 
Subroutine GNRATE passes the vector N to subroutine COUNT. 
(N is the 2* dimensional vector containing elements no, a, 


being the number of times, out of s trials, that z, zes, 


Seecurs.) 
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Let B denote an index of the elements in the vector N. 
Suppose .B = 4 represents both receivers 1 and 2 detecting 
a common signal, and over many signals this combination 
occurs twelve times. This event is then annotated. as 
N(4) = 12. Note that N(1) is the zero vector and cannot 
be observed. A convenient method of arranging B is to 
consider the binary number system. Figure 1 demonstrates 
how the binary numbers can be used to order B for a 
4-receiver system. For example, if receivers 1, 2, and 3 
all detect a common signal, this 1s associated with the 
bomeaay number 0 1 1 1 and the related value of B is the 
base ten equivalent, 7. B relates directly to the random 
vector Z. By computing a value for B, Z can be obtained 
simply by computing the binary equivalent of B. One 
convenient algorithm for computing an outcome on Z is as 
follows: 

1. Let B= 1 

2. Drew a U(0,1) random number 
Se ra U(0,1)sT,, fete se= 5. gic} 
4. If isk, go to 2 and repeat 


se Let Z = base.B, 


Z 


Using Figure 1 as a reference, it can easily be shown 
how the necessary counting is performed. Notice that the 
array Of binary numbers is actually a reflection of the way 
each receiver is listed by the index. By considering the 


k 


binary array to be a (2° X k) matrix, various row and column 


sums can be used as multipliers in the counting procedure. 
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N( B) Index Receiver i, i=1,4 Associated 
4 


1 2 3 Binary Number 
eG) 0 Oro GeO 
nr 2) 1 1 Om O wl 
N( 3) 2 Oar. “0 
N( 4) 3 1 OmOi =| 
N( 5) 4 5 Oe 
N( 6) 5 1 3 Oe Ca 
N( 7) 6 3 0110 
N( 8) 7 1 5 Ol 1 a 
N( 9) 8 4 0 0 0 
N (10) 9 1 4 i OmeOies| 
N(11) 10 4 ion 0 
N (12) in 1 4 ee 
N (13) ee 3 4 i 1 2ORO 
N (14) 13 1 5 4 flier Ona 
Neel 5) 14 2 8 4 le @ 
N(16) 15 1 3 4 peeled 


Figure 1* Binary Counter for a 4-Receiver System 
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The total number of signals detected, t, 1s the sum 
over elements Nyreee DL in the vector N. The number of 
signals detected by recelver i, m,, can be computed by 
looking in the appropriate column in the binary matrix. 
Receiver 1 is represented in the array by the far right 
column. The value for m4 1s computed by summing those 
elements in the vector N for which a 1 appears in the far 
right column. This procedure is repeated over all k 
columns. 

The value of m*, m* = ties (the number of signals de- 
tected by receivers i and j simultaneously), is obtained 
by summing those elements in N which contain two or more 
ones in the associated binary number. Remember in com- 
Plena m , that if, for example, three receivers detect a 
common signal, then that element of N will be counted three 
times. Likewise, ina four receiver detect, that element 
of N will be counted six times. Hence, depending on the 
number of ones in the associated binary number, there is 
an appropriate multiplier assigned to each element of N. 


The value for 25 25 (mym,) is the sum over the totals 


@eeained in each column, done for all ie<«j. 


PeeeCcOMPUTING ESTIMATES: "SUBROUTINE SOLVE" 
Subroutine COUNT passes all the necessary accounting 
information to subroutine SOLVE which performs the opera- 


tions and optimizations required to produce estimates of 


at 





TT; and s as discussed in section III. Using the input 
parameter METHOD in the main routine, it is possible to 
Weacify which one or more methods to be used. 

1. Method 1: Maximum Likelihood 

Maximum likelihood estimators will be found using 
numerical methods because as k increases, Lagrangian and 
other techniques quickly become too burdensome. In this 
simulation, the non-linear maximization approach used is 
the Generalized Reduced Gradient Technique (see Lasdon, 
(Ref. 4]). This computer method is available in most 
computer libraries and is easy to use. 

Library routines DATAIN, GRG, and OUTRES are called in 
sequence. Routine DATAIN reads all necessary data while 
routine GRG performs the optimization. It is necessary for 
the user to supply a subroutine GCOMP which evaluates con- 
straint and objective function values (see Appendix A). 
Routine OUTRES provides the user information on each iter- 
ation of the maximization. The quantity of output can be 
varied by a user inputted flag in DATAIN. In this simula- 
tion, the iterative data is totally suppressed. 

2. Method 2: Method of Moments 

This method is perhaps the only one of the three 
discussed which can be computed by hand. Estimates for 
TT. and s are simple functions of the data provided by 
subroutine COUNT and are computed as discussed in section 


ed, 
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3. Method 3: Least Squares 


Again, numerical techniques are necessary in 
solving least squares estimators. An initial value of 
P(A) = .1 begins the least squares minimization. The 
value .1 is chosen because of the potential presence of 
trivial solutions when the system intercept probability 
mals within the interval [0,.1) (see Knorr, 1979). P(A) 
is incremented in steps of .01 on the interval [.1,1]}. 


~~ 


For each value of P(A), P (j) is computed using (4) and 


Ps (j) is numerically evaluated using the algorithm on 
page 19. The value of P(A) which minimizes the least 
squares expression (5) becomes the system intercept prob- 


ability and is used to compute Thy and $. 


POUL PULT: “SUBROUTINE OUTPUT" 

Two forms of subroutine OUTPUT are used depending on 
which part of the simulation is being performed. Recall 
that. in the first part of the simulation, 72 environments 
were examined one time each to obtain a feel for the be- 
havior of the estimators. A sample listing of the results 
is shown at Figure 2. 

The second part of the simulation focuses on analyzing 
the mean absolute error (MAE), bias and variance of the 
estimators for TT, and s. Twelve of the environments were 
selected for analysis and replicated 100 times each. For 


each replication, the MAE of the element in the vector T 


Za 





was computed by comparing the estimated values to the 
expected values of mM. The expected value of 1, is easily 
Memputea given l., the specified ability of receiver i to 
detect a signal of strength V, and the distribution of V. 
Given that the form of the function for We = T, (V) = 

1. ° Vv, E( TT.) = E(1, Se = 1, ¢ E(V). The expected 
value of a Beta(«,]) random variable can be shown to be 
x/(1+%). Hence, E( TT, ) Sl, «/(]1+<). The MAE of the 
receiver intercept probability becomes ABS ( TT, - E(TT,)). 
The MAE's are summed across the entire probability vector 
to produce an MAE estimate for the system. The 100 system 
MAE estimates are then sectioned into ten blocks of ten 
replications each, making use of the computer library 
routine SECTN. In this manner, estimated sample parameters 
On the system MAE (i.e. mean, variance, etc.) are computed 
along with estimated parameters of the sample parameters. 
Figure 3 shows a typical output listing for subroutine 
SECTN. (Remember that these sample parameters are on the 
estimate of the MAE for TT, not on the estimate of A di- 
mectely. ) 

Following the sample parameters is a listing of the 
initial input to the signal generator. Next, the method of 
analysis used is shown, followed by direct estimates of the 
mean and variance for each element in ‘TY. 


The identical procedure is followed for analyzing the 


Eemimator S. A typical output listing for the analysis on 
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S is shown at Figure 4. The above procedure is duplicated 


for each of the three methods discussed in section III. 


\ 
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INPUT PARAMETERS ICG SIGNAL GENERATOR 
* NUMBER OF SIGNALS GENERATED: 1900 
NUMBER CE RECEIVERS: 3 
ALFA PARAMETER! 100.0 
RECEIVER PROBABILITY VECTOR: 0.20 0.50 0.80 
Bere CGUNTS 
N( LL) = 22 ; 
Nt 2) = 83 
hia yee ee 
a 4) = ale 
oie] Te 
Nf 3) = 324 
N 7) = BS 


GeeMIZATICN CCEFFICIENTS 

GBSERVEO NUYBER OF DETECTIONS BY RECIEVER 1: 
1$7.0 507.0 197.0 

OBSERVED NUMBER OF NO-DETECTS BY RECIEVER I: 
719.0 409.0 119.0 ° 

CBSERVED NUMBER OF TIMES J RECIEVERS DETECTED A COMMON SIGNAL: 
417.9 413.0 86.0 

NUMBER OF SIGNALS INTERCEPTED: 916.0 


CPTIMIZATICN METHOD: MAXIMUM LIKELIHOOD/CONDITIONAL MAXIMUM LIKELIHCOD 


ESTIMATES OF INTERCEPT PROBABILITIES: 
0.20 0.51 0.80 


ESTIMATE OF THE TOTAL NUMBER OF TRANSMISSIONS: 992.79 


APTIMIZATICN METHCD: METHCD OF MOMENTS 


ESTIMATES OF INTERCEPT PROBABILITIES: 
0.20 0.51 90.81 


ESTIMATE OF THE TOTAL NUMBER OF TRANSMISSTONS? 985.05 


OPTIMIZATICN METHCOs LEAST SQUARES 
are lcM INTERCEPT PREBABILITY: G92 


ESTIMATES CF INTERCEPT PROBABILITIES: 
Ce20 0.51 0.80 


ESTIMATE CF THE TOTAL NUMBER OF TRANSMISSIONS: 995.65 ‘ 


megure 2. Sample Output - Form 1 
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; Vo eoaNAGYoES OF RESULTS 


The purpose of this section is to demonstrate the 
results of the simulation described in section IV and to 
discuss, in general terms, trends and observations about 
the behavior of the different methods of analysis as they 
apply to the model. The first part of the simulation used 
each of the three methods in 72 different environments to 
obtain estimates of the total signal population size. [In 
the second part of the simulation, twelve environments were 
replicated 100 times each and mean absolute errors (MAE) 
of the estimates for radio intercept probabilities, T+ 
and total signal population size, s, were obtained. The 
mean and variance of the MAE estimates were computed over 
all replications. Additionally, the mean and variance of 
the estimators themselves were computed. 

The data shown in Table 1 and Table 2 summarize the 
results obtained in the first part of the simulation. [In 
Table 1 are the estimates of signals generated for a 3- 
receiver system using each of the three methods of analysis. 
The input values for the radio intercept capabilities and 
the alpha parameter are listed in the first column. The 
actual numbers of signals generated are listed in the second 
column comprising an environment. The 6-receiver intercept 


systems are shown in Table 2. In the 6-receiver environments, 


5 





two receivers share a common value for 1,, thus, a vector 
of .6/.7/.8 implies that two receivers were specified to 
have 1, values Gmmicemeewo Of ./,8and two of-.8 . 

This part of the simulation was designed to give a 
feel for how the estimators would behave over a wide 
range of conditions and to help select an interesting set 
of environments from which to analyze in more detail. Re- 
source limitations precluded a complete analysis of all 
environments. Those situations indicated with an (*) are 
the ones selected for replication. 

Though no formal claims can be made about the estimators 
in the first stage of the simulation, a review of the data 
in Tables 1 and 2 seems to indicate that all the methods 
perform better as alpha increases (as the receivers become 
less dependent) and as the receiver intercept probabilities 
are increased. This result is intuitively appealing 
because changes in the parameters which increase the ability 
of a Signal to be detected or a receiver to detect a signal 
will result in an increase in signals detected. Asa 
higher portion of the total signals are detected, more data 
are available for estimating the total number of signals 
present. 

The results of the second part of the simulation are 
shown in Tables 3 - 8. Tables 3 - 5 show mean and variance 


estimates for the receiver intercept probabilities for 


S65 





alpha equal to 1, 10, and 100 respectively. Tables 6 - 8 
show mean and variance estimates for the total signal popu- 
lation sizes. Notice that all mean estimates of receiver 
intercept probabilities tend to be positively biased. The 
bias is more significant for alpha equal to 1 and improves 
greatly for alpha equal to 10, and tends to become nega- 
tively biased for alpha equal to 100. Thus, it appears 
that a high degree of dependence between receivers will 
have a significant impact on the bias of the estimators. 
Similar trends are not as readily apparent on mean esti- 
mates of population size. Though maximum likelihood and 
method of moment estimators seem to be negatively biased, 
the least squares estimator shows several instances of 
strong positive bias, particularly for small sample siZe. 
One possible explanation for this may be found in the large 
variance estimates for the least squares estimators. This 
may be an indication that the model is not adequately dis- 
tinguishing between the trivial and non-trivial solutions 
for P(A). There appears to be consistently much better 
estimates for signal population size as alpha increases. 
With the wide variety of changing parameters and rela- 
tively few replicated trials, it is not possible to make 
comparative claims between the estimators under changing 
conditions, but only to look at the estimators under the 
conditions tested. However, because there were four 


environments tested for each value of alpha, there are some 
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suggestions about the quality of each estimator which can 

be made. For each value of alpha, the sum of the means 

and meriances of the MAE estimates was computed over all 
environments tested for each estimator. These sums repre- 
sent a score for each estimator over various conditions but 
with constant alpha. These mean absolute error scores for 
estimates of radio intercept probabilities and total signal 
population size are shown in Table 9. Notice that the 
method of moments estimator appears to perform the best for 
estimates of when alpha equals 10 and 100. The variances 
of the estimates are, however, too large to claim it is the 
best estimator. This table does show that there is a sig- 
nificant reduction in error as receiver dependence decreases. 
No one estimator appears to be better for signal population 
Size over all tested values of alpha, but again, the esti- 
mators do tend to improve as alpha increases. It is inter- 
esting to note that maximum likelihood estimates tend to 
have the smallest variability while the least squares esti- 
mates have the largest. 

Ideally, all 72 environments should be replicated and 
estimates compared over all conditions to obtain performance 
criteria on each estimator. A data analysis of this mag- 
nitude would provide valuable insight into the problem but 


is beyond the scope of this thesis. 
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Input Signals Signals Maximum Method of Least 
Values Generated Detected Likelihood Moments Squares 
x,l, Estimates Estimates Estimates 
= | 
2S 19 19.40 18.74 19.59 
may. 7/.8 100 ee I Si 5 Gore 3 TonO0 
1000 ie 805.01 TA Ga GS S358 
10000 65 8248.49 7510509 8384.21 
ZS 14 14.98 14.62 15256 
me/.3/.4 100 49 Gili 17 53270 62e05 
1000 566 747.70 698.84 Tad 
10000 5614 7886.43 7512.48 7907.04 
ZS 17 AN) 012 Zine 137554 
27 .5/.8 100 70 150 | 74.83 To Z 
1000 719 126205 7a eso 138 
10000* 7344 8104.96 MG06. 7 2 8251.68 
oc = 10 
25 ZS 25.04 24.87 Zou 0 
my 277.8 100 100 1COnm 23 99.54 100.00 
1000 9916 Oe 993 eke 996.200 
10000 99169 O99 Z 8 9955.48 996 300 
25 24 Zo 2733 280,97 
aay. 3/.4 BOG. G1 104.85 LOS oo 105.81 
1000 844 Sito: | 978.60 981.40 
10000 8568 GID 2 9866.30 9962279 
Z5 24 24333 24.44 Zoe 
0/8 100 97 98.40 oor cS O33 
10101 990 1006.67 gg SS 2 1010320 
10000 oe oo Dove 9924.80 9910.10 
ec = 100 
25 oS Zae0 | 24.85 2004 
hoy 27/.8 TOO oD SS) Oy oon 9900) 
1000 998 9398.94 Soe? epeeits) 8 (C)0) 
10000 9994 TOCO0R SS 9992.41 9994.00 
z= Ze Ze a2 Zoe 6 26 os 
eyes 7. 4 100 79 SGe¢ 94.81 86.81 
1000 876 985.43 978.01 9954 
10000 * 8777 9307 7.50 Jeger os 9973.86 
25 24 2442 Zar. 93 24.00 
2 ao . 8 100 oo oe uOG7.7 0 101.02 
1000 999 1006.50 O07 247 999.00 
10000 9926 NEOO2 23 D976. 55 9926.00 
Table 2 ° Estimates of Signals Generated/6-Receiver System 


(Two receivers share a common intercept probability) 
40 





l 


ALFA 


100 


NUMBER OF STGNALS GENERATED 


3 


NUMBER OF RECEIVERS: 


1009 


NUMBER OF SIGNALS GENERATED 


3 


NUMBER OF RECEIVERS: 


co 


NUMBER OF SIGNALS GENERATED 


6 


NUMBER OF RECEIVERS 


1832 1.882 


DAODUUADA OO UW 
MOAOM SPL 
eee eos eos” 


41 


10000 


NUMBER OF SIGNALS GENERATED 


6 


NUMBER OF RECEIVERS: 


OOWNOOU 
aa St 
e©eee © oS 


VER INTERCEPT PROBABILITIES. 
OVER THE ENTIRE 


UTED 


ECEI 
COMP 


TABLE 3. 





10 


ALFA 


22 


NUMBER OF SIGNALS GENERATED: 


3 


NUMBER OF RECEIVERS: 


©0441 1.223 


10000 


NUMBER OF SIGNALS GENERATED: 


3 


NUMBER OF RECEIVERS: 


00 


e 
e 


NUMBER OF SIGNALS GENERATED 


6 


NUMBER OF RECEIVERS: 


TaN DOW! 
Am OO ama 
Madr oom 
aI AIMCO CY 


cd 


| aD ff OO LL 
att VALCO Ll 
eee eees 


Se 


42 


10090 


NUMBER OF SIGNALS GENERATED 


6 


® 4 
e 


RECEIVERS 


NUMBER OF 


BABILITIES. 
IRE 


TABLE 4. 





109 


ALFA 


25 


NUMBER OF SIGNALS GENERATED: 


5 


NUMBER OF RECEIVERS: 


1000 


NUMBER OF SIGNALS GENERATED: 


3 


NUMBER OF RECEIVERS: 


100 


NUMBER OF STGNALS GENERATED: 


6 


RECEIVERS: 


NUMBER OF 


rm Ot Olt ~F Os 
Ot Oath 
DX HODOAA 
LAVIN OO RMN 


on 


OD BROMUS 
WIA OO P< 
eeeecs of 

bol 


43 


10000 


NUMBER OF SIGNALS GENERATED: 


6 


NUMBER OF RECEIVERS: 


Og OO 
NALCO CN PP 


oeeee oS 
—— 


OBABILITIES. 
TIRE 


PR 
EN 


PL! 


TABLE 5. 





l 


ALFA 


M Lh 
— LU 
NY <f 
It = 
LLJ = = 
ao 


ENERATED 


27 


8 
le 


(MAE) 


+t WO 
DOO 
e «E> «@ 
OOO 
To © 

oOo mom 

wf 
OOOO 
fLHoed 
6 @ ¢ 8 
Om 
mo wy 
rf NN 
OTFO 
OM ta 
ee @ 6 
uy wy 
NO + 

. 

LLJ 

) 

me 

<— 
a 
Ieaszx< 
LL) <{ ~~ <{ 
2S 

© 

NS 

e 

~ 

LY 

— 

e@ 

~ 

) 
(T) md 

i 
\t 

= 
ice 


on 
min 


ro 


OO 
hr 


UVOI 

e °@ 
mt 
Oo 


ANCE 


ae 
Iw 
LL <f 
=> 


(MAE) 


1000 


N 


a tO 
ONIN = 

oS e e@.hU6©¢ 
OmO00 
mt OLA 
NMA 


UNO OLY 
oncr 
eet e@ 
OO ANNO 
le @2@ soe) 
NPN 


IND COW 
Taio 
eeee 
CONWAY 
ND OW 
NOAA 


ANCE 


ae 

Ie 
LL <f 1 <f 
= oz 


3 
©30/.35/ «40 


I 


K 
p 


44 


OO 
al) =] 
ee 


Ow 


m+ 
mil 
e 6 


(MAE) 


uy moO 
LAN UN © 
eeee 
NU UN 
Tm oO 

CO a= 


TOAW 
DaDo 


iF © 


Om tC 
Door 
eee *@ 


mr © 


ANCE 


as 

Irae 
LLJ <{m~ <f 
=o = 


Q0/.35/.35/ 


tr 
o~ 
aS 
cot 


Or 
OO 


~o 


ANCE 


ae 
ax 
Lu <f 
=> 


| 


(MAE) 


Hoe 
= . oe Pec acd/ 


Na ov) 
eoeee 
CM WyO0 
—M 00 
M~F OO 
taf 


ON MAN 
© @e @ 
NATE oO 
Ono aor 
LV at (} SO 
NMOAIN 


Om Os 
e@eee 
DOF I 
~ONM 
Ont CO — 
a yest 


ANCE 


| 

qv Zé 
LL <f =~ <f 
= 


10/.10 
40/.40 


out 


I 


N 
K 
p 





10 


ALFA 


S) LJ 
fue LL) fe 
Nat 
I= 
LL — 
=Jd J 


GENERATED 


43.42 


3136 


OO — 
OM 
e °® 


ro 


MEAN 
VARIANCE 


(MAE) 


20.28 
LO=> 
205.00 


3083.2 


an © 
DOUICO 


NPFOO 
=~ 


ODOM 
FUVE OS 
e © @ @ 


NOY © 


ANCE 


a 

{Irw~znx< 
LL) <f{ + <f 
oe. 


°18/.45/ .73 


mt 
—= CO 
Ow 
OD —4 


LA 
mt Lt) 
OO —4 
Qi = 


Win 
ro 


(MAE) 


10000 


=aHOdm 
OPN 4 
® @s6 @ 
mat PUY 
oo wu 

Oo rw 


U\ tend 
SODO 
e@ @ @ ¢@ 
ttor 


Our wo 
“+ xt 


Our Ph CO 
rt ON P- CO 
ee @ e 
UV et i0\ <O 
“Ph AN 


ANCE 


a= 

Ita 
LJ <f 1 <f 
zo 


3 
©0554 264/ 213 


45 


i -—~ 


© @ 
eae e) 
NN 


or 
pe net 
se @ 


Ou’ 


FO 
uVOS 
e 
OS = 
rm 


ANCE 


<c 
LL <f 
=> 


(MAE) 


POW 
Omi oO 
* eO e 
PNA 

ot ee 


O*CO-e1 
u—VO eO 
ee e 
LVUN OLN 
am i) 


OS <P 
FO 
e* @ @ ®@ 
SY | 
~~ md 


$e2U/.20/ 


996.8 
42.66 


e 
sO 
ru 


mo 
» @ 
(7 
Onl 


VARIANCE 


MEAN 


(MAE) 


AUN pO 
Om | 
ee *@® ¢ 
LAV emt OD 
an | 


set OO 
FoOolf 
e © e 
DOO 
+r MO 


KXFUO 
NOcal 
oe 0 @ 


FN WY 


ANCE 


ae 

Ie ax 
LL) <f{ 1 <f 
> > 


1645/.45/ 


8 
3 


ION SIZE. 
NG THE 
ENERATED.) 


bm oe 9) 





100 


ALFA 


bh LU 


ENERATED 


ANCE 


=— = 
I 


= 


(MAE) 


5 


= we 
tro;o 
e © » 
CQNaW 
Na 
= 


Om © 
NaNO 


ONOW 


6.69 
16.98 
46.86 


ANCE 


ae 

ava 
LL <f mw <t 
Lo T= 


~20/.30/ .40 


3 


“AO 
or 
ON 


ANCE 


ae 
Ic 
LL <X 
=> 


(MAE) 


1000 


N 
K 
p 


DAFA 
very 1 N 
ee «© e 
NI Oan 
aM + 


FAD 
me OOO 

on en 
QM 4 
AMO W 


ANCE 


a 

Iemzx 
US <{ = <f{ 
oe ee aa 


3 
020/.50/.79 


| 


46 


100.6 
wee 


100.00 
© 93 


ANCE 


ae 
Iw 
LL} <x 
=> 


(MAE) 


COMO 
Or 1 © 
e eC e@ 

el (T) 


ONO tr 
COU es 
ee «© e 


OND 
“OO 


ANCE 


a 

qIeawx 
Us <f mt <f 
=> T= 


Oe 


Ow 
Ov 


MEAN 
VARIANCE 


(MAE) 


uy MO 
IFNaO 
ee ¢@ @ 6 
OM 
OMNI AY 
— — 


aa UN 
FAIA 
oe ee 
COS AN 
FUN 

GJ = 


mes 
CO LN CO 
ee @ © 


AON 


N 
TANCE 


<Icazx< 
LL <f m1 <f 
27-2 


eal al 


02 
04 


aA. oO 


See 


Citi 
MWS 





MEAN ABSOLUTE ERROR SCORES 
ON ESTIMATES OF RADIO INTERCEPT PROBABILITIES 


Maximum Method of Least 
Likelihood Moments Squares 

Estimates Estimates Estimates 
mona = 1 2. E (MAE) 1.5272 1.8873 1.7616 
> Var (MAE) .0611 .0993 sie 5 
Alpha = 10 2ZE(MAE) ~5904 -5684 .6677 
2Var (MAE) B02 9 O23 7 .0996 
Alpha = 100 2 E(MAE) ~6046 15822 5905 
2 Var (MAE) .0261 O29 5 .0449 


MEAN ABSOLUTE ERROR SCORES 
ON ESTIMATES OF TOTAL SIGNAL POPULATION SIZE 


Maximum Method of Least 
Likelihood Moments Squares 
Estimates Estimates Eseimates 


Alpha = 1 SXE(MAE) 2241.54 2787.95 2115.60 
X= Var (MAE) 4591.70 4465.26 14003.57 
Alpha = 10 XE(MAE) 66.34 102.47 121.84 
=Var (MAE) 598.71 Wigs 4107.67 
Alpha = 100 SE(MAE) 117.85 68.81 94.11 
Var (MAE) 2461.50 1437.97 2915.05 


Table G* Mean Absolute Error Scores 
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VI. SUMMARY AND RECOMMENDATIONS 


This thesis has presented a development of various 
estimators which can be used to estimate the total number 
of signals passing through an acquisition system given 
only a portion were detected. Likewise, estimators were 
found for estimating the intercept probabilities for each 
receiver in the system. These estimators were used in the 
context of a model developed to represent the signal de- 
tection process. The model suggested is the foundation 
for a simulation conducted over a variety of conditions 
including receiver and signal characteristics as well as 
transmitter activity. The simulation results provided some 
general insight into the quality of the estimators devel- 
oped and their behavior over varying conditions. The 
question of dependence between receivers was briefly ad- 
dressed and it was shown that as dependency increased, the 
quality of the estimates decreased. 

The work presented here is not viewed as a final solution 
to the problem. The model used was deSigned for flexibility 
and it would be a simple matter to refine and expand the 
detection process used within this report. The simulation 
conducted was sufficient for validating the model and ob- 
taining some initial insight into the quality of the esti- 


mators but was not adequate for making conclusive claims 
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about them. It is hoped that this thesis may stimulate 
others into continuing work on this subject. In particu- 
lar, it 1S suggested that a complete experiment be designed 
and a formal data analysis of the simulation results be 
conducted to determine which method of estimation is opti- 
mal. It would equally be interesting if someone with an 
electronics background could incorporate into the model 
additional parameters or functional relationships bringing 


in some hardware considerations or radiowave characteristics. 
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APPENDIX A 


SIMULATION PACKAGE CHARACTERISTICS AND REQUIREMENTS 


The simulation discussed in section IV is designed to 
Minimize user interaction. The user must decide what is 
desired to be accomplished by the simulation and modify 
the main routine and output routine accordingly. The 
material discussed in this appendix suggests simple in- 
structions on how to use the FORTRAN simulation package 
provided with this study. A complete source listing of 


the FORTRAN code is located at the end of this appendix. 


mee =CHARACTERISTICS 

As currently designed, the simulation is dimensioned 
to handle an acquisition system of up to ten receivers. 
The total number of signals which can be generated in any 
one replication is 99,999 , limited only by the input for- 
mat statement allowing five digits. Similarly, the random 
number generator seed is formatted for up to five digits. 
To modify the simulation to handle a larger acquisition 
system or increase the number of signals desired to be 
generated is a simple matter of increasing the dimension- 


ality of the arrays and enlarging input formats. 


50 





eee Mo lHOD SELECTION 

The*user must determine which one or more methods of 
estimation to be performed. Should only one method of 
estimation be desired, that method is requested through a 
coded input parameter in the data base as follows: 


Maximum Likelihood....METHOD 1 


Method of Moments..... METHOD 2 


eas ErogUanes ss... Mi LHOD on 

If more than one method is desired, the user should provide 
the lowest coded input parameter in the data base and 
establish the appropriate looping in the main routine. The 
loop established in the source listing at the end of this 


appendix providesfor all three methods to be applied to a 


common signal data set created by subroutine GNRATE. 


Ce REQUIREMENTS WITH THE USE OF GRG 

If the user specifies the maximum likelihood method of 
estimation to be used, the user supplied subroutine GCOMP 
must be included in the source listing. This subroutine 
contains the objective function to which the generalized 
reduced gradient (GRG) routine attempts to maximize. (In 
actuality, the GRG routine minimizes an objective function, 
thus, to maximize the log-likelihood function it is neces- 
sary to minimize the negative of the objective function.) 
Subroutine GCOMP contained in the source listing provided 


is general in nature and accurately represents the log- 
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likelihood function discussed in section III for any 
k-receiver system. No user modifications to this routine 
are necessary for systems of ten receivers or less. 

Each time the maximum likelihood method is used during 
the course of a simulation, the GRG routine is called. 
With each call to subroutine GRG, it 1S necessary to pro- 
vide a special input data base to support the optimization. 
This data base is read by library routine DATAIN, and is 
supplementary to the initial data base established for the 


Simulation. 


D. DATA INPUT 

Six data cards initialize the simulation. An addition- 
al set of data cards are necessary for each call to sub- 
routine GRG. The size of the data set required for GRG 
1s related to the number of receivers. Shown at Figure 5 
is a sample data set for a three receiver system. The 
data is preceeded by a data section number and followed by 
the input format. Below, is a description of each data 


entry by data section number. 


Data Sctn Number Data Elements and Description 
1 Random number generator seed (read once). 
2 Number of signals to be generated. 
3 Number of receivers in the system. 
4 Input receiver intercept probability 
vector. 


SZ 





10 


vel 


12 


3 


14 


eS 


16 


yy 


Value of abpha used in generating B(« ,1) 
variates. 


Method of analysis to be used. 


Problem name. Any combination of 80 
alphanumeric characters. 


Number of receivers (values of 1. to 
be estimated) followed by a '1' anid 
'O', the number of equality and in- 
equality constraints in the maximiza- 
EAlrOn. 


Number of variables with finite lower 
bounds (equal to the number of receivers). 


There should be k cards each containing 
a pair of numbers. The first number is 
the index of the variable in 9 above 
followed by its lower bound. [In this 
Simulation the lower bound is set at 
.0005 


Number of variables with finite upper 
bounds (again, equal to the number of 
receivers). 


There should be k cards each containing 
a pair of numbers. The first number is 
the index of the variable in 11 above 
followed by its upper bound. [In this 
Simulation the upper bound is set at 
~9999. 


Number of inequality constraints having 
finite upper bounds. Enter 0. 


inilelat- values for TY, Bnitwer- e700 for 
each. 


Tolerance value when certain iterative 
processes should stop. Default value 
is zero. 


Limits on the upper bound of allowed 
iterations. Default value is zero. 


Specification on the desired amount of 
output. Zero is the default value and 
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a value of one allows the user to in- 
crease or decrease output. A value of 
one is used to suppress output in this 
Simulation. 


18 A pair of zero values effects the sup- 
pression of output. 


Je Anything but 'QUAD' will assure the use 
of tangent vector extrapolation for est- 
imating initial values of basic variables, 

20 If the user does not wish to partially 
or completely specify an initial basis, 
enter a zero on this card. 

For a more detailed discussion on GRG refer to a GRG user's 


guide (i.e. see Lasdon, 1975). 


Mee OUTPUT SELECTION 

Two forms of subroutine OUTPUT are provided with the 
source listing. The first form is used when one or more 
environments are specified and replicated once each. This 
form would typically be used with actual field data. The 
second form of subroutine OUTPUT is used when it is desired 
to replicate one environment a number of times to allow 
for a statistical analysis of the estimators T, and 6. 
The user need only to decide which form of the output is 
applicable to the simulation being conducted and insert 


the appropriate listing. 
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Data 
Section 
Number Sample Data Format 
Card Column 
ipaooon oo idle i) io 2222222223333333333 
012345678901234567890123456789 


i PAA 5 £5 

Z 10000 es 

5 3 Le 

4 "eo, .70 .80 MOS ec 

5 100). ES | 

6 1 15 

i Problem name. Any 80 Alphanumerics. 20A4 

8 5 1 0 BES 

9 3 ES 

10 1 70005 Too ol Ov. 4 
2 .0005 
3 .0005 

a 5 5 

2 1 ~9999 155 57 5 1064 
2 ooo 
3 ogee 

13 0 ns 

14 “a0 220 ~oC 8E10.4 

1. 0 nS 

16 0 Tes 

Vy 1 ILS 

18 0 0 215 

19 0 A4 

20 0 nD 


Figure 5 ¢ Sample Data Set for a 3-Receiver System 
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APPENDIX B 


LIST OF SYMBOLS 
(alphabetically) 


A: the event (Z#0) 
B: an index of the elements in the vector N 
tJ: number of data vectors z for which Zz. = j 


feet?! /(n,!...n ) 
: ok_4 


k: number of receivers 
L(s,p): likelihood of N at n 


1.: proportion of signals of strength V that receiver 1 
would detect 


M;: number of signals detected by receiver i (outcome is 
m., ) 
ms 2-; ., (number of signals detected simultaneously by 


receivers 1 and j) 


N: random vector with outcome n; NaM te (S70) 


2 °=1 
Mee(T) , Nlaypsse, Ni me Vveccor Of n.*S in some specific 
O 1 k Z 
2-1 order 
n,: number of times z S occurs 
Mem, D,;<-e, DP ): vector of p_'s in same order as 
O 1 k Z 
2-1 fOr mh 
P, oj): mass function of z-. 


Stes 2; /'1 - P,) 
r=2 ° - 1 


S: sample space of 2 
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S$: number of signals 


S,7, 3 .estimates of s and ,, respectively 


fey.) = MW... = Ll. © V. 
ay 1j z 5 


t: number of signals detected 
U: uniform(0,1) random number 
wo 2 


strength of signal j (outcome is v0 


J 
Z= (Z.,,..-, 4): vector of Bernoulli detection or non- 

1 k ee . 
detection indicators, for a given signal, 
for all receivers 

Zz = a : umb f j hi 
_ = j=1 243 Number of receivers w ich detect a given 
Signal 
re. probability receiver 1 detects signal j 
ee S(T. (V)) = 1, * ECV) 


«: parameter of Ve where ae 
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